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Abstract. Based on additional interferometric observations, we reanalysed the CO(2-l) and 3mm continuum 
emission of Q0957+561, a lensed QSO at a redshift of z = 1.4141. The emission in the CO(2-l) lines reveals a 
gas-rich host galaxy with a peculiar double-peaked profile at one of the two lensed images. Our new interferometric 
CO maps of the host galaxy agree well with HST images obtained by Keeton et al. (2000) and we thus argue that 
the two velocity components arise from molecular gas in the disk of the host galaxy. We also present new model 
calculations, all in excellent agreement with recent time delay measurements and simulations. 
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1. Introduction 

Since the discovery of Q0957+561, the first confirmed 
gravitationally lensed quasar (Walsh et al. I1979f) at a 
redshift of z = 1.4141, several models have been devel- 
oped to understand the lensing potential of the interve- 
ning galaxies: a giant elliptical galaxy (Gl) at a redshift 
of z = 0.355 with a surrounding cluster and probably an- 
other group of background galaxies at z = 0.5 (Angonin- 
Willaime et al. [TUMI Chartas et al. ITMSjl . Besides flux 
amplification, gravitational lensing helps also to improve 
spatial resolution as detailed by Kneib et al. (1998). 
Differential lensing effects which are detectable through 
the comparison of lensed images at different wavelengths 
indeed provide a powerful tool to probe the structure of a 
galaxy at a much higher resolution than possible with cur- 
rent millimeter interferometers. Planesas et al. (1999, in 



the following P99) have succeeded in detecting the CO(2- 
1) line in the host galaxy around Q0957+561, but lacked 
a detailed lensing model to investigate the source of mole- 
cular emission. Keeton et al. (2000) have recently detected 
also the stellar component of the host galaxy showing dif- 
ferences in its distribution with respect to those of the gas 
obtained by P99. 

To confirm and substantiate the weak line profiles pre- 
viously detected by P99 and to eliminate the previous in- 
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consistencies to Keeton et al. (2000), we have carried out 
new observations with the IRAM interferometer. To fur- 
ther improve on P99's interpretation, we have developed a 
numerical code which incorporates existing lensing models 
of Q0957+561 1 . 

2. Observations 

CO(2-l) and (5-4) observations of Q0957+561 were made 
simultaneously in 1998 and again in 2003 with the IRAM 
interferometer. The first observation run is discussed in 
P99. We only give additional details on the refined re- 
duction of this data set. 

Data taken in 1998: The observations were carried out 
in the B, C and D configurations. The bandwidth of the 
cross-correlator was set up to cover 420 MHz at both fre- 
quencies with a spectral resolution of 2.5 MHz, equiva- 
lent to 7.8 km s- 1 at 95.5 GHz (redshifted CO(2-l)) and 
S.lkms" 1 at 238.7GHz (CO(5-4)). The bandpass was 
calibrated on 3C273, 3C345 or 0923+392 while phase 
and amplitude calibrations were performed on 0917+624. 
Unlike P99, we discarded data from an entire observing 
run that were substantially reducing the quality of the 
maps. We have also corrected for the velocity scale which 
was wrong by a factor of 1.5 in P99 due to a data header 
error. 

Data taken in 2003: The observations were carried out 
in the C and D configurations. The bandwidth of the 



1 We adopt a flat cosmology with f2M=0.3, Qy=0-7 and a 
Hubble constant of 65 kms~ 1 Mpc _1 . 
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Fig. 1. Continuum at 3mm (left panel), integrated CO(2-l) line emission (middle panel) and the spectral CO(2-l) 
profiles observed toward CO— A, CO— B and at the 3mm position of the radio jet (right panel). Offset position (0,0) is 
at the assumed position of the lens, i.e. at aj2ooo =10:01:21.0, <Jj2000 =+55:53:50.5. Contour levels are from 5<r = 1.1 
to 4.7mJybeam _1 in steps of 0.6mJybeam _1 for the 3.1mm continuum and from 3er = 0.33 to 0.77 Jykms -1 in steps 
of 0.11 Jykms -1 for the line emission. The synthesized beam of 3.2"x3.0" at p. a. 87° is shown in the lower left corner 
(1" corresponds to 9kpc at an angular-size distance Da — 1.87 Gpc). The white stars indicate the optical positions 
of the A and B image. The rms noise per 10 MHz channel is ~0.65mJy. The spectrum at the position of the jet is 
indicative of the noise level. Velocities are relative to the redshift at z = 1.4141. 



cross-correlator was set up to cover 580 MHz with a spec- 
tral resolution of 1.25 MHz corresponding to 3.9 km s -1 
at 95.5GHz and 1.6kms -1 at 238.7GHz. We resampled 
the two data sets to a spectral resolution of 7.8 km s" 1 
(2.5 MHz) over the 420 MHz band covered by the first ob- 
servations to ensure uniform noise in the merged velocity 
channels. 3C273, 0851+202 and 0420-014 were used as 
bandpass calibrators while 0804+499 and 1044+719 were 
taken as phase and amplitude calibrators. Relative to P99, 
the sensitivity of the combined data sets is now higher by 
a factor of ~2. 

3. The data 

3.1. Continuum emission 

3.1mm: The radio continuum was computed separately 
for the two data sets by averaging the line-free channels 
with a velocity lower than — 510 km s -1 and higher than 
+540 km s -1 . Two lensed images of the quasar, labeled 
A and B, and a radio jet C appear at this wavelength 
(FigQJ, in agreement with VLA observations (Harvanek 
et al. I1997f) . The positions coincide with the optical ones 
within the errors (Tab. 1). Also in agreement with previous 
work, the A and B components are pointlike whereas the 
radio jet C component is slightly elongated (~ 2.7", PA 
^60°) in the beam of the interferometer. Finally, there is 
no evidence for variability above 10% between May 1998 
and April 2003 in all three components. This coincides 
with results obtained at 3.6cm with the VLA by Harvanek 
et al. l(TM7|> . 



1.3mm: A noise level of a = 0.6 mJy was obtained in the 
2003 data by averaging over the entire bandwidth available 
at 1.3mm. We find a 3.5<r peak at the 3mm position of the 
A-continuum component, a 2a peak at the position of the 
B component and a 4<r peak at the position of C. The 
steep drop in flux density at 1.3 mm is consistent with 
the weakness of the synchrotron emission expected if we 
extrapolate to this short wavelength the spectral index of 
~ —0.6 measured in the 20 mm (Harvanek et al 1997) to 
3 mm (P99, and this paper) range. Any dust emission is 
either too weak or too extended (say > 3") to be detected 
at this frequency. Data from the 1998 run were not taken 
into consideration as they were much less sensitive. 

3.2. Line emission 

CO(2-l): To estimate flux densities in the CO(2-l) line 
(Tab. 1), the 3.1mm continuum was first removed sep- 
arately from each data set to account for possible low 
level continuum variability. Then, all data sets were com- 
bined to a single uv-table, and velocity channels in the 
— 510 km s -1 to +540 km s -1 range was summed to pro- 
duce an integrated line map (Fig^). The resulting map 
shows two lensed images labeled CO— A and CO— B (P99). 
The emission centroids of the two images are separated by 
~ 7", about 1" more than in the optical, radio and mil- 
limeter continuum (Tab. 1). The spectral profiles taken 
towards positions CO-A and CO-B are different. While a 
double-peaked profile is visible towards the northern im- 
age CO-A, a single blueshifted velocity component is de- 



M. Krips et al.: Q0957+561 revised: CO emission from a disk at z = 1.4 



3 



CO(2-»l) Continuum at 3.1mm 6 

A-blue A-red B-blue ABC 



Rel. positions (") 


(-1.9,43.8) 


(-1.6,4-6.1) 


(40.2,-1.1) 


(-1.45,45.00) 


(-0.24,-1.0) 


(3.23,7.13) 


Position errors (") a 


(+0.2,4-0.2) 


(40.3,40.3) 


(40.2,40.2) 


(40.06,40.08) 


(40.08,40.10) 


(40.07,40.08) 


Flux density (mjy) 








4.740.2 


3.4 4 0.2 


6.2 4 0.2 


Peak flux c (mjy/beam) 


0.9±0.2 


2.140.2 


2.240.2 


4.640.2 


3.340.2 


5.140.2 




(1.1±0.3) 


(1.740.3) 


(2.340.3) 


(4.240.3) 


(2.540.3) 


(5.540.3) 


FWHM (kms- 1 ) 
Vel. at line peak d (kms -1 ) 


280±60 


160420 


280450 








-140+30 


340415 


-140423 








Integr. line flux (Jykms" 1 ) 


0.25±0.06 


0.3440.06 


0.6140.06 









Table 1. Flux densities, velocities and positions of individual image components. Values have been determined by 
fitting Gaussian profiles to the visibilities. Offset position (0,0) is at the assumed lens position specified in the caption 
of Fig^ a Relative position uncertainties are based on statistical noise and do not include astrometric errors. 6 Values 
taken from the second observation campaign in 2003. c Values in brackets are from the reanalysis of the P99 data, and 
do not include data from 2003. d Velocities are relative to the redshift at z = 1.4141. The CO(2-l) line parameters for 
the two A components were not derived from the spectrum shown in Fig^lthat was taken between A-red and A-blue, 
but are based on spectra determined at the respective positions of A-blue and A-red. 



tectcd towards CO-B. The shape of the blue components 
measured towards CO-A and CO-B is similar. The CO— F 
image reported by P99 could not be confirmed with the 
new observations. 

CO(5-4): No significant line emission was detected at 
the position of CO— A and CO— B, the two line compo- 
nents detected at 3mm. A ~ 1.9 mjy (= 3er) peak is 
only tentatively detected in the velocity integrated map 
(-300 km s" 1 to +100 km s" 1 ) at the position of CO-B. 
The intensity of the CO— B component at 3mm over the 
same velocity range is estimated to be 1.8 mjy. Although 
uncertain because of the low signal to noise ratio of the (5- 
4) line, and possibly because of some residual continuum 
emission, we set an upper limit of < 1 on the velocity 
averaged intensity ratio R54 =CO(5-4)/CO(2-l) toward 
CO— B, the strongest of the two line components. 

4. Modelling Q0957+561 

We developed a numerical code based on the standard 
gravitational lens equation (Schneider et al. I1992fl to ex- 
plain the absence of the double-peaked line profile toward 
CO— B and to investigate the distribution and kinematics 
of the molecular gas around Q09574561. 

The code has been applied on 6 different models 
of Q0957+561 (Tab. 2), all based on previous work by 
Barkana et al. 1199j| and Chae et al. (1999). The mass 
distributions of three models are based on King profiles, 
as proposed by Falco et al. iftUgSl - models: FGS, FGSE, 
FGSE+CL), the remaining three on a softened power-law 
distribution, as suggested by Grogin et al. (1996-models: 
SPLS, SPEMD, SPEMD+CL). FGS and SPLS are spheri- 
cal models where the effects of the surrounding cluster at 
z = 0.355 are approximated by an external shear. FGSE 
and SPEMD models take into account the ellipticity in the 
lens galaxy, the remaining models use a single isothermal 
sphere (SIS) to model the lens properties of the cluster 
(SPEMD+CL, FGSE+CL). A point mass was added to the 
King profiles to account for a black hole in the center 



of the lens galaxy (Mediavilla et al. 120001 Barkana et al. 
1999), except for the softened power law models which im- 
plicitely cover this case. The group of background galaxies 
at z = 0.5 has not been considered, however. The compos- 
ite pseudo-Jaffe models by Keeton et al. (2000) have not 
been taken into account as they would not have provided 
further details for our analysis. 

To find the set of parameters for each model reproduc- 
ing our observations with the lowest Xr > we nrs t restricted 
the parameter space with simple assumptions based on 
the number of lensed images, separations among them etc. 
We then scanned the parameter space for each model to 
get the lowest \r and included the following constraints: 
relative positions of the A and B continuum components 
at 3.1mm with respect to the lens (4 constraints), rela- 
tive positions of CO- A-blue and CO-B-blue (TabEJ, the 
continuum and line intensity ratios and the optical/radio 
time delay of 400-420 days from Kundic et al. (1997) and 
Haarsma et al. (1999). This results in a total of 11 con- 
straints. For simplicity, compact fixed-size Gaussians were 
used to approximate the respective components. The so 
found best-fit models turned out to be in excellent agree- 
ment with simulations based on VLA and optical data 
(e.g. Barkana et al. 1999, Chae et al. 1999 and Keeton et 
al. 2000). 

The SPEMD+CL, FGSE+CL and the FGS models re- 
produce the observed constraints with the lowest Xr (< 3; 
Tab. |2J). Although the FGS model yields one of the best 
results, we do not favour spherical mass distributions be- 
cause of the ellipticity of the lens (Bernstein et al. 0.997). 
The contribution of the cluster is also important: The 
models that best explain these observations all require an 
SIS cluster (FGSE+CL and SPEMD+CL, the latter has a 
lower x 2 ) 2 ■ We note that the time delay derived from the 



2 For the following discussion, we will mainly use the best- 
fit FGSE+CL model since it provides the closest match to 
already published results, e.g. Keeton et al. (2000). However, 
both models are fully equivalent for argumentation purposes. 
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SPLS e SPEMD e SPEMD+CL 1 FGS^ FGSE e FGSE+CL e 



Source Positions 
(Source plane) 

Radio continuum 

(Aa s ,AS s ) (") a 
Line emission 

(Aas,A5 s ) (") a 


(-0.02,0.90) 
(-0.23,0.95) 


(-0.14,1.02) 
(-0.32,1.00) 


(2.20,2.49) 
(1.93,2.39) 


(-0.15,1.06) 
(-0.38,1.05) 


(-0.18,1.10) 
(-0.40,1.05) 


(1.80,2.58) 
(1.54,2.47) 


Lens properties 
Galaxy: 
a v (kms -1 ) 

e 

$PA (°) 6 

V 

Cluster. 

i d 

%> (°) " 
er c i (kms 1 ) 

(AQcl.AdcO (") 6 


400 
0.12 

1.1 

0.3 

60 

— 


400 

0.12 

0.1 

64 

1.0 

0.2 

60 

— 


360 
0.02 
0.2 
64 
1.1 

375 
(13.7,6.9) 


300 
1.2 

0.4 
60 

— 


340 
0.6 
0.1 

64 

0.3 

60 

— 


320 
1.1 
0.3 

64 

- 

350 
(13.7,6.9) 


Black hole: 
M bh (1O 12 M ) 








0.5 


0.3 


0.2 


Ar cont [days] 
Arii nc [days] 9 

_ total 

■m RC 


411 
421 

5 

6 


420 
450 

4 

7 


417 
412 

4 

9 


417 
423 

3 

5 


421 
454 

3 

6 


416 
414 

3 

6 


iVdof c 


3 

3.0 (3.4) 


2 

5.5 (40) 


2 

1.5 (1.5) 


3 

1.3 (1.2) 


2 

5 (40) 


2 

2.0 (2.5) 



Table 2. Best fit parameters for Q0975+561. Parameters are: Act = right ascension, AS = declination, a = velocity 
dispersion, £ c = core radius, e = ellipticity, <I>pa = position angle, r\ = power law index, 7^ > = shear term, $ 7 = shear 
angle, Mbh = mass of the black hole, At = time delay, m = magnification factor. a Offset from assumed lens position. 
b Fixed parameters (shown in italic). c Values in brackets are Xr - tests with time delay. The reduced x? is defined as: 
Xr = X 2 /A r dof where iVdof is the degrees of freedom. d 7' is only taken for the elliptical potentials (7' = 7/(1 — k)) 
where k = S c i us /S cr i t is defined as the surface mass density of the cluster in units of the critical surface density 
(Barkana et al. 1999). For the spherical models, k = and for the elliptical ones it can be approximated by k ~ 7, i.e. 
7' = 7/(1 — 7). e based on Barkana et al. (1999). ^ based on Chae et al. (1999). 9 determined for the blue component. 



same models are in rather good agreement with recent op- 
tical and radio measurements (~ 400 — 420 days; Kundic 
et al. 1997; Haarsma et al. 1999). 

4.1. Line emission 

In the following two subsections we discuss different ap- 
proaches for modelling the blucshiftcd line emission based 
on the best-fit FGSE+CL potential and then explain the 
difference seen in the profiles toward CO— A and CO— B. 

4.1.1. The blueshifted line 

Let us first concentrate on the blueshifted line. We can 
propose two different scenarios for the origin of blueshifted 
line emission: 

— a compact or extended 3 region centered on the nucleus, 

3 Simulated by an elliptical Gaussian profile, i.e. a distribu- 
tion which is symmetric with respect to its centroid. 



— a compact or extended region spatially not coincident 
with the nucleus. 

The small line ratio A-blue/B-blue of ~ 0.4 already rules 
out a compact region centered on the nucleus, since it dis- 
agrees with the measured 3mm continuum ratio A/B of 
1.5 ± 0.1. Even if we account for a contribution of the jet 
in the A component and assume a lower ratio for the core 
emission at 3mm, we know from VLBI observations that 
the ratio A/B at the core is close to 1.3 (Garret et al. 
1994). Thus, the continuum at 3mm and the blueshifted 
CO emission cannot originate from the same compact re- 
gion (compare also Fig|2J). 

On the other hand, an extended line emitting region 
centered on the nucleus can indeed result in an A-blue/B- 
blue of 0.4. To constrain the size of the blueshifted region, 
we assumed an elliptical Gaussian profile and fitted the full 
widths at half maximum (FWHM) and the position angle 
(PA) depending on the continuum and line ratios between 
the lensed images. A best-fit (Fig|3Jl is obtained for a disk- 
like region, that is extended in the direction of the nuclear 
jet (~ 5") and tiny in the direction perpendicular to the 
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CO: 



— 2 - CO: different positions for blueshifted emission 



-2.5 - 




RA (offset) 



Fig. 2. x 2 results for the continuum (upper panel) and 
blueshifted line (lower panel) position fits based on the 
FGSE+CL model. Shown here are different positions in the 
source plane for the blueshifted line emission and the con- 
tinuum depending on the % 2 -test including the observed 
intensity ratios between the two lensed images and the ob- 
served position. Assumed values for the 3mm continuum: 
A/B= 1.5; blueshifted line: A/B= 0.4. Positions for the 
A, B (continuum) and CO— A and CO— B can be found 
in TabQ] Contour levels are from 0.1 to 1.0 in steps of 0.1 
(in terms of exp(-x 2 ))- The region with a value close to 
1.0 in this map corresponds to the best-fit position. 




-FWHM ["] 



0: blueshifted component at off — nuclear position 



2 - 




1 2 3 

x-FVVHM ["] 

CO: red + blueshifted component at nuclear position 



jet (~ 0.2", PA=20°), with the south-western part of the 
disk crossing the inner tangential caustic of the lens. A 
priori, such a model produces multiple images (up to 5; 
compare also Keeton et al. 2000) that result in a stronger 
CO— B image. 

Alternatively to the extended region centered on the 
nucleus, blueshifted emission from a rather compact region 
(< 1") close to the inner tangential caustic and slightly 
off (~ 0.5" to the south-west) the nuclear position is also 
consistent with a line ratio of ~ 0.4 (Fig|2J). A symmetric 
extended region (> 1") at an off-nuclear position does not 
produce such a low A-blue/B-blue of 0.4. 

Thus, an independent blueshifted emission must either 
originate in an extended region centered on the nucleus or 
in a compact region with a position close to the tangential 
caustic. 

4.1.2. The redshifted line 

The difference in the line profiles observed toward CO— A 
and CO— B can be explained by the location of the red- 
shifted gas component relative to the lens caustic. The 
best-fit simulation for the FGSE+CL model produces one 
tangential caustic (FigQJ . The radial caustic is at infinity. 
While the continuum and blueshifted line components re- 
main close to the tangential caustic and are thus deflected 
into two distinct images, the redshifted component is al- 




V'HM ["] 



Fig. 3. x 2 results of the size fit for the line emission. The 
FWHM of the Gaussian profile are plotted with the re- 
spective x 2 - We also varied the position angle (PA) from 
5 to 60° but we restricted the presentation to the PA with 
the lowest x 2 • Upper panel: Blueshifted emission centered 
on the nucleus, x 2 values based on the A-blue/B-blue ra- 
tio. Lowest x 2 for PA=20°. Contour levels start at 0.56 
to 1.0 in steps of 0.04. Middle panel: Blueshifted emission 
centered at a fixed off-nuclear position closer to the tan- 
gential caustic ( Exposition of line emission in Tab|2I). x 2 
based on the A-blue/B-blue ratio. Lowest x 2 for PA=20°. 
Contour levels are from 0.1 to 1.0 by 0.1 (in terms of exp(- 
X 2 ))- Lower panel: Red- and blueshifted emission centered 
on the nucleus, x 2 values based on the (A-blue+ A-red)/B- 
blue ratio. PAs of 50±10° all result in a low x 2 - Contour 
levels are from 0.6 to 1.0 in steps of 0.08. 
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Right Ascension (offset) ["] 
5 0-5 




Fig. 4. The line-free continuum at 95 GHz (upper left panel) and the 12 CO(2— >1) line emission (upper middle panel) as 
observed with the IRAM interferometer from both data sets. The lower panels show the simulated radio continuum (left) 
and line emission (middle) convolved with the synthesized beam and the critical lines (grey lines). In the right panel 
the location of the unlensed line emissions (r = A— red,b = B— blue = A— blue) as well as the unlensed radio continuum 
(J = jet, Q = quasar) is shown together with the radial and tangential caustics with respect to the assumed position of 
the lensing galaxy (grey lines). The offset position (0,0) of each map corresponds to the position of the lensing galaxy. 
Contour levels for the data are as Fig^ We have also used the same contour levels for the simulated maps. The black 
stars indicate the optical positions of the A and B image. The parameters used for the model (FGSE+CL) are listed in 
Tab. [5] Changes in the velocity dispersion of the lens of ~ 1% or in the source positions of ~ 0.1" produce a difference 
in the time delay of the FGSE+CL model of — 5%. 



ready located too far off the caustic (to the north-east) to 
generate two lensed images. Thus, it appears that both 
blue- and redshifted gas become visible toward CO— A 
while blueshifted gas only can be detected toward CO— B. 
The unlensed separation between the centroids of the blue- 
and the redshifted components, corrected for magnifica- 
tion, are estimated to or equivalently ~ 9kpc at the 
distance of Q0957+561. 

Likewise, an extended region centered on the nucleus 
where the north-eastern part corresponds to the redshifted 
component and the south-western part to the blueshifted 
one produces equivalent results (FigEJ. To obtain the de- 
rived line ratio (A-blue+ A-red)/B-Blue ~ 1.0, the size of 
the region has to be extended by about ~ 3 — 4". This 
value seems to be independent of the lens potential even 
if we find that the SPEMD+CL model appears to favour 
slightly more elliptical profiles for the emission than the 



FGSE+CL model (cf. also Keeton et al. 2000). If we di- 
vide the lensed CO— A into a red- and blueshifted part 
and derive then the ratio between A-blue and B-blue, we 
again obtain a value of ~ 0.4 4 . 

The sensitivity is not sufficiently high to establish the 
absolute genuineness of the CO— A-bluc component from 
the observations alone, but our best-fit FGSE+CL model 
of the host galaxy not only agrees with the Barkana model 
of Keeton et al. (2000) , but a priori corroborates the detec- 



4 In contrast to Section l4.1.1l where we found that an ex- 
tended blueshifted emission region would not result in an A- 
blue/B-blue ratio of 0.4 if symmetrically centered on a position 
closer to the tangential caustic, the extended emission region 
for the blueshifted part is not symmetric with respect to the 
position used there. This might explain why we nevertheless 
obtain such a low ratio for the blueshifted part alone. 
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Fig. 5. HST image in the H-band of the host galaxy 
in Q0957+561 (Keeton et al. 2000) with contour levels 
from Fig0]of the velocity integrated CO(2— 1) emission. 
Separation of the tick marks is 1". 

tion of blueshifted CO gas in the direction of component 
A. 

4.2. Radio continuum 

The morphology and the intensity ratios of the radio 
continuum can easily be modelled. As already noted by 
Schneider et al. I|1992[) . the radio jet J is lensed only once 
due to its location far off the caustic. 

5. Discussion 

The 3mm continuum observations of the redshifted dou- 
ble quasar Q0957+561 are in good agreement with obser- 
vations at other wavelengths. While the 3mm continuum 
emission from the quasar is unresolved, emission from the 
radio lobes is extended along the direction of the jet. The 
CO(2-l) line, however, does not originate in the quasar. 
The difference in position, and therefore in the deflection, 
explains also the difference between the flux density ratios 
in the 3mm continuum (A/B = 1.5 i.e. like in the radio 
continuum) and in the line (CO— A/CO— B = 0.4 for the 
blue component), and accounts for different time delays 
for the CO(2-l) and the radio images. We note, however, 
that the time delay difference is only marginally significant 
(at most 5—20 light days, according to the models). 

What can we conclude about the origin of the red- 
shifted CO(2-l) velocity component? P99 argued that 
it was either originating from a companion galaxy, very 
close to Q0957+561 (see also Papadopoulos et al. [2001) 
or tracing the presence of molecular gas in the rotating 



disk of the host galaxy. Our models discard the first 
hypothesis. First, there has been no detection of such 
a companion galaxy at other wavelengths, but the host 
galaxy, which has already been revealed in HST observa- 
tions by Keeton et al. (2 0"0"0"|l . is oriented in the direction 
of the two line components in CO— A and appears to be 
extended by ~ 3" (Fig. Furthermore, models for the 
gas emission in the host galaxy that take the velocity 
profiles observed toward CO— A and CO— B into account 
are similar in size and shape (see Fig0J to those invoked 
by Keeton et al. (2000). The blue- and redshifted CO(2-l) 
emission appear to be very likely connected to the stellar 
distribution of the host galaxy and thus corroborate the 
results of Keeton et al. (2000), whereas models based 
on a single blueshifted component favour a very thin, 
elongated disk or a compact CO (2-1) emitting region 
slightly off the nuclear region, and as such are rather in 
contrast with images of the host galaxy. Moreover, the 
integrated line ratio (CO-A-red+CO-A-blue)/CO-B 
~ 1.0 is consistent with the ratio H-band- A/H-band-B~l 
derived by Keeton et al. (2000). This is what one would 
expect if the profile of the gas distribution is very similar 
to the stellar one. Also, the double-peaked line profile 
at CO— A which appears to be symmetric and centered 
within the errors at z = 1.4141 is a classical signature for 
rotation and therefore, provides further support that the 
CO emission is associated with the host galaxy. By the 
same arguments we are tempted to exclude also a recent 
galaxy merger in Q0957+561, though this possibility 
cannot be ruled out completely. 

Based on these arguments, we favour P99's second 
hypothesis: the presence of an important reservoir of 
molecular gas in the disk of the host galaxy surrounding 
Q0957+561. 

According to the models the continuum is less magni- 
fied by a factor of 2—3 relative to the strongly lensed opti- 
cal images (~ 10) while the total magnification factor de- 
rived in the CO line varies between m^Q 1 ~ 6 (FGSE+CL) 
and m£?o al — 8 (SPEMD+CL). Based on the lower magnifi- 
cation, we have estimated an upper limit to the molecular 
gas mass for the integrated blue- and redshifted veloc- 
ity profile in CO-A and CO-B. The FGSE+CL model 
corresponds to the higher CO— luminosity case and pro- 
vides therefore a slightly higher molecular gas mass. Under 
the assumption that the brightness temperature of the 
(1-0) and (2-1) lines is the same, we adopt a CO to H2 
conversion factor of ~ 4.8 M (Kkms _1 pc 2 ) -1 based on 
a determination for the Milky Way (Solomon & Barrett 
[TMTIl . Clearly, the Milky Way and the Q0957+561 host 
galaxy are in different evolutionary stages and show differ- 
ent properties, so the adopted conversion factor is likely 
to be an overestimate (cf. Downes et al. 1993), and an 
upper limit for the molecular gas mass will be obtained. 
Based on the FGSE+CL magnification factor, we estimate 
M gas =M(H 2 +He)~ 2.4 x 10 10 M Q for the blue profile and 
- 2.7 x 10 10 M Q for the reddened profile (~ 9 x 10 9 M Q 
for SPEMD+CL), giving an upper limit for the total gas 
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CO-A-red CO-A-blue CO-B-blue 



mL co 








(10 10 Kkms~ 1 pc 2 ) 


0.8±0.3 


0.8±0.2 


1.8±0.2 


wfgsb+cl b 


1.5 


1.7 


4.3 


M gaa (10 10 M s ) 


2.7±1.0 


2.4±0.6 


2.3±0.6 



total gas mass Mg£ = MgH + M^ c ~ 5 x 10 10 M, 



Table 3. Apparent CO luminosity L' and molecular gas 
mass M gas of Q0957+561. a m = magnification factor. 
b determined with the FGSE+CL model. 



mass of ~ 5.1 x 10 10 Mq. We can also give a lower limit 
for the gas mass of M gas ~8x 10 9 M Q assuming that the 
CO(l-O) line is optically thin (Solomon et al. HW7|l . This 
relatively small difference of only a factor of 6 indicates 
that the standard conversion factor might nevertheless al- 
ready give a reliable estimate of the gas mass in the host 
galaxy around Q0957+561. Assuming a radius of ^lOkpc 
and a velocity separation between both line profiles of 
~ 400km s -1 , we find a dynamical mass of ~ 4 x 10 11 M Q 
(not including inclination effects). This is about ten times 
higher than the derived gas mass consistent with what 
is found for other active galaxies (e.g. Evans et al. 2002, 
Sakamoto et al. 1999). The low upper limit on the velocity- 
averaged line intensity ratio, ~ 1, favours the hypoth- 
esis of low excitation CO emitted mainly from the disk of 
the host galaxy. The agreement between molecular gas 
masses obtained with the integrated CO luminosities and 
individual magnification factors, each tracing line emis- 
sion from half of the quasar host, is further support for 
the rotating disk hypothesis. If our assumptions are cor- 
rect, the molecular gas in the disk cannot be very hot, but 
to our knowledge no sensitive observations in the ground 
and higher rotational CO transitions have been carried 
out yet to confirm this conclusion. 

6. Conclusions 

Recent sensitive observations were combined with earlier 
data by Planesas et al (1999) to corroborate their original 
discovery of CO(2-l) line emission from Q0957+561. A nu- 
merical program was developed to analyse the properties 
of the lensed system, the results of which can be sum- 
marized as follows. While the region of blueshifted line 
emission is found to lie in between the two caustics, and is 
therefore deflected into two images, redshifted emission is 
found outside the caustics, and therefore results in a single 
lensed image. We further argue that both redshifted and 
blueshifted line emission originate from the same system: 
a disk with a molecular gas mass of ~ 1 — 5 x 1O 1O M in 
the host galaxy in Q0957+561. Our results for the host 
galaxy are in excellent agreement with a previous optical 
work by Keeton et al. (2000). To our knowledge HR10 
(Andreani et al. I200U|I and Q0957+561 are the only sys- 
tems at redshifts of 1 < z < 2 in which CO emission was 
clearly detected, but in contrast to Q0957+561 where we 
have not yet been able to detect dust emission, HR10 is 



rich in dust and molecular gas. Although both systems 
are significantly different, they both independently give 
crucial insights into an epoch during which the bulk of 
stars of the present day Universe formed, and thus mark 
an important phase in galaxy evolution. 
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